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Abstract 

The relation between Loop Quantum Gravity and Regge calculus has been pointed out many 
times in the literature. In particular the large spin asymptotics of the Barrett-Crane vertex 
amplitude is known to be related to the Regge action. In this paper we study a semiclassical 
regime of Loop Quantum Gravity and show that it admits an effective description in terms of 
perturbative area- Regge-calculus. The regime of interest is identified by a class of states given 
by superpositions of four-valent spin networks, peaked on large spins. 

As a probe of the dynamics in this regime, we compute explicitly two- and three-area corre- 
lation functions at the vertex amplitude level. We find that they match with the ones computed 
perturbatively in area- Regge-calculus with a single 4-simplex, once a specific perturbative ac- 
tion and measure have been chosen in the Regge-calculus path integral. Correlations of other 
geometric operators and the existence of this regime for other models for the dynamics are 
briefly discussed. 

PACS: 04.60. Pp; 04.60. Nc 

1 Introduction 

The relation between Loop Quantum Gravity [H [21 [3] and Regge calculus [Ij has been pointed out 
many times in the literature, both at the kinematical and at the dynamical level. 

The fact that in SU (2) representation theory some inequalities and symmetries appear which 
have a geometrical interpretation in terms of flat n-simplices is known since long ago [6l [7| . Its 
relevance here comes from the key role played by SU (2) representation theory in Loop Quantum 
Gravity [S] on the one side, and by flat n-simplices in the Regge approximation to General Relativity 
[3] on the other. For instance, restricting attention to kinematics, four-valent nodes of a S'C/(2)-spin- 
network share some geometrical properties with flat tetrahedra and can be considered their quantum 
version [lOl [TTl [12] . Moreover, appropriate superpositions of these 'quantum tetrahedron'-states 
can be shown to describe exactly the geometry of a classical tetrahedron [T3] . 
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Spin foam models |14U15l[TB] provide a covariant way to code the dynamics within the framework 
of Loop Quantum Gravity. Despite the fact that the spin foam formahsm can be motivated from 
various independent perspectives, building a model always requires some input. For instance, this 
input can come under the form of a classical action principle and a procedure to make sense of its 
formal path integral. In any case, once a model is available, one can forget the heuristics and start 
studying its properties. In order for a specific model to be claimed to be a model for quantum 
gravity (and not to be any other diffeomorphism invariant theory) some check of its semiclassical 
properties is necessary. 

The most well-studied model for quantum gravity is certainly the one proposed by Barrett 
and Crane in [T7j. It is given by a sum over intertwiners of the {15j} symbol squared. More 
recently Engle, Pereira and Rovelli [IHl [TH] proposed an improvement from the Loop Quantum 
Gravity point of view of the Barrett-Crane model. Both models take as starting point some form 
of classical simplicial gravity and end up with an elementary "symbol" in SU (2) representation 
theory. However this by itself is not enough to guarantee that the dynamics admits Regge gravity 
as a semiclassical regime. An indication that in fact this is the case, at least for the Barrett-Crane 
model, was presented in [20] where it was shown that the Barrett-Crane vertex amplitude has a 
large spin asymptotics of the form 

^Bc(jmn) ^ P{jmn) COS {Sncggc{jmn) + ^ ) + D{jmn) (1) 

where S'Roggo is the Regge action for a single 4-simplex with faces having areas related to the spins 
jmn- This result is analogous to Ponzano and Regge asymptotic formula for the Wigner {6j} 
symbol giving the cosine of Regge action in three dimensions . It can be considered as a formal 
semiclassical limit, with h —^ and the areas kept fixed. 

Expression ^ was obtained starting from an integral formula for the Barrett-Crane vertex 
|21j and using a stationary phase approximation. In [20j . Barrett and Williams shown that the 
stationarity condition for the phase reproduces exactly the form of the Schlafli differential identity 
satisfied by the Regge action (i.e. the discrete version of Palatini identity). This identifies the Regge 
action up to an overall scale. Among the stationary points there are both ones corresponding to 
four dimensional non-degenerate simplicial geometries and ones corresponding to low-dimensional 
degenerate and non-degenerate geometries. The latter ones contribute to the D{jmn) term in ([1]). 
Numerical computations [22] and analytical calculations [23l [24] show that in fact the degenerate 
configurations give the dominant contribution to the asymptotic formula. This has raised doubts 
on the physical viability of the model. 

However it should be noticed that simply taking the large spin limit is not enough to identify 
a semiclassical regime, exactly as taking the large distance limit in the Hydrogen atom transition 
amplitude kernel (r2|L''(T)|ri) is not enough to identify its semiclassical regime. What is needed 
is a state peaked at large distance from the attractive center and on a specific momentum. The 
state codes the initial conditions for a classical orbit and the dynamics of this state identifies a Ke- 
plerian semiclassical regime as can be shown computing correlations of position at different times. 
Taking into account this fact, in [25l[26] a study of the dynamics of the Barrett-Crane model on a 
boundary semiclassical state peaked on a specific intrinsic and extrinsic 3-geometry was presented 
and the following general mechanism was shown to be at work: of the right-hand side of ([T]) only 
the exponential of i times the Regge action contributes to the two-spin correlation function, while 
the other terms are suppressed by the phase of the semiclassical state. In this paper we strengthen 
and extend this result. 
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The main initial motivation for the present work was to compute the three-spin correlation 
function at leading order in the framework described above. A task of this kind requires some 
difficulties to be overcome the first. In [251 I26j . in order to compute the spin-spin correlation 
function needed for the analysis of the graviton propagator, formula ^ was widely used. However 
the detailed dependence of D{j„m) and P{jmn) on the fluctuation of a large spin jo was not given 
in [231 [21] and it cannot be deduced from numerical computations [22] ■ In US US] D{jmn) was 
assumed to be slowly varying with the fluctuation Sjmn = jmn — jo and P{jmn) was assumed to 
be independent from the fluctuation Sjmn, at least at leading order in j^^. Moreover, the same 
function 5'Rcgge(jmn) is not completely under control as it is given as a function of the spins which 
are related to the area of the faces of the 4-simplex [571 [251 [IH] and not to its edge lengths: besides 
the fact that 

- for a single 4-simplex there are isolated configurations for which giving the areas does not 
completely fix the geometry of the 4-simplex, and 

- for more than one 4-simplex some constraints between the areas arise which are still to be 
worked out, 

the main problem is that an explicit expression for the deficit angles in terms of the areas is still 
not available. Nevertheless, the relation between area-Regge-calculus and Loop Quantum Gravity 
pointed out by Rovelli in [30] seems to be both inevitable and desirable, as acknowledged in [5T] . 
The class of difficulties discussed above can be completely avoided if one manages to compute the 
spin-spin correlation function directly from the Barrett-Crane vertex formula, without using its 
asymptotic expansion llj. This was done in [32] where a suited boundary semiclassical state was 
introduced in order to perform analytically the calculatior0. The result confirms the scaling found in 
[25l[26] and puts it on stronger grounds. However the analysis in [32] presents two disadvantages: (i) 
it hides completely the relation with area-Regge-calculus at any intermediate step by construction 
and (ii) the new semiclassical state is not suited for a phenomenological parametrization beyond 
the leading order. 

The route we follow here is to separate the two problems. As in |32j . we study correlation 
functions in Loop Quantum Gravity without using expression ([!]), but still using the semiclassical 
state used in [251 126| which allows a phenomenological parametrization for corrections contributing 
to the three-spin correlation function. This is possible thanks to the two following results: 

• writing the Barrett-Crane vertex as an integral over the group manifold (5*^)^, we identify 
the region of the manifold which gives the dominant contribution to the n-spin correlation 
functions. We find that it is given by a bounded domain B of radius of order The 
integral on its complement gives an exponentially suppressed contribution; 

• the dominant contribution can be computed explicitly thanks to the properties of the two 
parameters coming into play, the large mean spin jo and the small relative fluctuation fc„j„ = 
Sjmn/ joi as coded by the boundary semiclassical state. Hence we have both 

- an asymptotic expansion in jo^^ and 

- a perturbative expansion in A:,„„ 

of the Barrett-Crane integral restricted to the domain B. 

These results allow to compute again, with a different technique, the two-spin correlation function 
computed in |25[ I26[ [32] . Moreover, using this technique, we compute the three-spin correlation 

^For a similar analysis in three dimensions, see |33l 134) . 
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function. Its relevance for the study of the semiclassical regime of Loop Quantum Gravity will be 
discussed in a separate paper |35| . Here it is enough to notice that, as the two-spin correlation 
function is related to the graviton propagator, the three-spin correlation function is related to 
the 3-graviton Feynman-interaction-vertex. But there is much more: the analysis discussed above 
identifies a perturbative regime within the non-perturbative Loop Quantum Gravity framework. 
Thanks to the boundary semiclassical state, the dynamics can be described perturbatively in the 
fluctuation Sjmn in terms of a perturbative action and a perturbative measure. 

On the other hand we study perturbative quantum area-Regge-calculus and compute two- and 
three-area correlation functions for triangles belonging to a single 4-simplex. By perturbative we 
mean that we fix a background Regge skeleton and background values of the areas so that they 
describe flat space. Then we quantize perturbatively the fluctuations of the areas in a path-integral 
fashion. The quantum theory is defined in terms of a perturbative action and a perturbative 
measure. This approach was developed in [36l [37] for length- Regge-calculus. Starting with a single 
regular 4-siniplex as background, area variables pose no problem as the Jacobian from lengths to 
areas is well-defined. While the perturbative action can be obtained from the standard Regge 
action, the measure in the integral over area-fluctuations has to be postulated. Many proposals 
have appeared in the literature for the non-perturbative measure over lengths and, although the 
specific measure does not affect area correlation functions at leading order, it does at higher orders. 

This approach of keeping Loop Quantum Gravity with semiclassical states on the one side 
and perturbative quantum Regge calculus on the other allows a direct comparison between the 
two frameworks. We find that a complete matching can be achieved once a specific non-trivial 
perturbative measure has been chosen on the Regge calculus side. 

The paper is written in a self-contained form and is organized as follows. Sections [2][6] are 
dedicated to the Loop Quantum Gravity calculation. In particular in section [2] we introduce the 
framework and discuss our assumptions; in section [3] we give a representation of the Barrett-Crane 
vertex as an integral over (5*^)^; in section [4] we identify the region of {S^)^ in the integral which 
gives the dominant contribution to correlations and in section [5] we compute such contribution 
using an asymptotic expansion. We conclude the Loop Quantum Gravity analysis in section [6] 
where we give an effective description of the Barrett-Crane model in terms of a perturbative action 
and a perturbative measure. Moreover, in section [HI we give explicit expressions of two-area and 
three-area correlation functions on the chosen semiclassical state. Section [7] is dedicated to the 
perturbative quantum Regge calculus calculations. We set the framework for the calculation of 
area correlation functions and give an explicit perturbative expression for the area-Regge-calculus 
action. The comparison between the Loop Quantum Gravity calculation and the Regge calculus 
one is presented in section [7?2l which concludes our analysis. The last section of the paper, section 
\8\ is dedicated to a discussion of the result, to an analysis of some perspectives it opens and to the 
possibility of using the techniques introduced here to test the recently proposed new models for the 
dynamics [M [H [351 IMl SOI SI] • 

2 Large-scale correlations in spinfoam models 

In this section we introduce the setting for the calculation of correlation functions in Loop Quantum 
Gravity. The boundary amplitude formalism is widely used both here in the Loop Quantum Gravity 
calculations and in section [7] in the Regge calculus calculation. 

In the boundary amplitude formalism one considers regions TZ of the 4-manifold A4 and assigns 
(i) a Hilbert space of states Tis to the boundary manifold S = dTZ, and (ii) a map Wu : Hs C 
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which codifies the dynamics in the region TZ. This approach has been advocated in various forms in 
[1], in [42l|43l|44l|45], in [46] and in [471 HH |49] . This approach presents both technical and concep- 
tual advantages. On the technical side, it avoids the need to prescribe the semiclassical asymptotic 
behaviour at spatial infinity [50] . In a sense, it could be the route to properly define asymptotic 
flatness in quantum gravity. On the conceptual side, it shifts the attention from "quantum cos- 
mology" [51] to quantum gravity and the hope is that it could open the possibility to calculate 
scattering amplitudes in background independent approaches |52) . 

2.1 The boundary amplitude formalism at work 

Here we shall take a minimalist approach (the reader can refer to [IB] for the general philosophy). 
We consider 4-dimcnsional Riemannian Loop Quantum Gravity with the dynamics implemented 
using the spinfoam formalism. In particular we study the Barrett-Crane spinfoam model jl7j and 
some of its modifications [53], (TB]. We restrict attention to a subspace of the boundary Hilbert 
space, the space of states having graphs containing only 4-valent nodes. The Barrett-Crane model 
is well defined on this subspace. 

As well-known IJ, in Loop Quantum Gravity the Hamiltonian constraint acts on a spin- network 
state non-trivially only at nodes. In spinfoam models, the action of the Hamiltonian constraint at a 
node is given in terms of a spinfoam vertex amplitude [54j . In particular, the Hamiltonian constraint 
corresponding to the Barrett-Crane spinfoam model is supposed to act in the following way on a 
state with a 4-valent node 




where on the right-hand side we have depicted the spinfoam corresponding to it. Such transition 
amplitude can be written in the following way: 

( X 1^1 X ) = {E^fij^))Mjl,^n) ■ (2) 
/ 

To capture the role of the vertex amplitude, this formula is best written in the boundary amplitude 
formalism introduced above. The strategy is the following: (i) cut out from the two-complex a 
4-ball B4 containing a spinfoam vertex v as shown in the following picture 



5 



(ii) introduce a vertex amplitude Wy to codify the dynamics in the region B4; it is a map from 
the boundary Hilbert space Tigs to C; (iii) introduce a state ^"53 ,j[s] belonging to the boundary 
Hilbert space 7^53 to describe the state on = dB^. The role of the state ^53 ,j[s] is to codify the 
dynamics outside -B4. As the boundary of the 4-ball intersects the two-complex giving the complete 
graph with five nodes Um and ten links Imn, 




the boundary Hilbert space is in fact an Ti.r^ which has the spin networks \ji2, ■ ■ ,i45,*i, ■ ■ ,^5) as 
a basis. Hence, instead of equation now we have 



W^Omnjin) = {Wy\ jJ^ifM ) = ( 11 ^/ 0™) ) 0™" ' ) ■ (4) 




Prescribing the action of the Hamiltonian constraint in this way, i.e. specifying the vertex ampli- 
tude, makes it easier to guarantee its crossing symmetry [54] . 

The quantities we are interested in here are correlations of geometric operator^ at the vertex 
amplitude level, such as area-area correlations on a state ^'rs.gOmn, *n) 



or volume- volume correlation {Vn' Vn" )q- Besides studying two-point correlations, we analyse n- 
point correlations too. Technically, these are correlations of coloring^ as for instance [55 l [56 l [57 l [58] 



Am'n''^r5,q{jmn,in) = SttG \/ jm'n' [jm' n' + l)*r5,g(jmn, in) ■ (6) 



^In the following, by Vn we mean the volume operator VCJZn) for a region TZn which contains the node n. Similarly 
we denote by Amn the area operator A{Smn) for a surface Smn which is cut by the link from the node m to the 
node n. These are well-defined operators on Ti-r^- 
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The motivation for computing correlations at the vertex amplitude level comes from the fohowing 
key fact: in order to have large scale correlations in a realistic situation where an appropriate 
boundary semiclassical state is chosen (and a sum over two-complexes is considered), correlations 
have to be present already at the level of vertex amplitudes. The absence of such correlations 
would provide an obstruction to the existence of a semiclassical regime. This fact was stressed by 
Smolin in [SS], where the ultra- locality issue was first addressed; see also [SD], [52 and [53] for 
a discussion. 

2.2 Spinfoam vertex amplitudes for quantum gravity 

Now we analyse some specific models for the dynamics in the Loop Quantum Gravity framework. 
They are defined giving a vertex amplitude Wyljmm'in), i-e- a map from Tips to C. The most 
elementary "crossing symmetric" vertex amplitude is certainly given by ^i;(jmn,*n) — {15j}sc/(2)- 
This vertex amplitude defines a diffeomorphism invariant quantum theory without local degrees of 
freedom. It corresponds to classical SU{2)-BF theory. 

Here we are interested in 4-dimensional riemannian gravity, which is locally 50(4) invariant. 
The idea then is to start with a S'0(4)-invariant weight for the vertex in the two-complex and to 
take its projection on a 5*0(3) subgroup of 50(4) on the boundary, i.e. on the graph labelling 
representations of 50(4) in terms of SU{2)l (8) SU{2)ii representations we have 

A wide class of models can be obtained using the most elementary 50(4)-invariant two-complex 
weight, namely the topological one given by a 50(4)-invariant BF theory 

where As!7(2) (jmn, *n) is given by a {15j} symbol. To guarantee crossing symmetry, the branching 
function /''^"''r r r can be assumed to be factorized over nodes in the following way 

rn<n n 

Both the Barrett-Crane model [17] and the Engle-Pereira-Rovelli model [18^ fall in this class. 

In the following we use standard diagrammatic methods for SU{2) representation theorj|f|. The 
{I5j}su{2) symbols in formula ([8|) can be written introducing a virtual link in representation k to 
decompose the 4-valent intertwiners in terms of two 3-valent ones: 

-4s(7(2)(jmn,fcn) = '/aI!" 'IJH ' ^^^^ 



^For an introduction see for instance appendix A of the book the lectures |64| or the book [65 1 . 
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Moreover, we take the following normalization for 3-valent intertwiners: 




1 if ii,j2, is are admissible, 
otherwise. 



2.2.1 The Barrett-Crane model 

We restrict attention to the Barrett-Crane model. This model is defined taking a branching function 
of the form ([9]) and defined in the following way: on links we have 

"^J'l." ^ (-1)''"" ^l,„r^,U^r^+3S^J 

where we have used the name Imn for the representations on the links of the boundary spin-network. 
On nodes (using the decomposition introduced in equation pU)) ) we have 

,/c„ _/ (-l)^''"(2fc^ + l)4,^,fc« for fc„ admissible, 
" [0 otherwise. 

Due to the 5jL jr and the 5^.L f.R^ in this model only simple 5*0(4) representations (on links, 
as well as on nodes) contribute to the vertex amplitud^. From p2|) one finds that the model is 
defined only on a subspace of Tirs, '^^^ with integer representations, i.e. it is well defined only 
for 5*0(3) spin-network states. In the following we call such representations Imm instead of jmn- 
The sign factor (—!)'""> will be justified lateJl. Due to the form of /^^ ^.h, the model is independent 

from the specific splitting of SU{2)l and SU{2)j^ intertwiners used in pO)) . Moreover, as f^l is 
independent of fc„, we have that this model weights all the admissible intertwiners with the same 
weight. 

The spinfoam boundary amplitude kernel defined in equation (|4]) can then be written as 

W4Un,in) = ( n (2/™n + l)'^^J(-l)^'""ABc(^mn/2), (14) 

with the face amplitude generally taken with exponent Nj = 2. The function A^c in (O has been 
introduced in order to make contact with the nomenclature generally used in the literature. We 
call it the Barrett-Crane vertex function. It was introduced in Tf] and is defined by the following 
formula 

5 

^Bc(jmn) ^ JJ(-l)^*='(2A:i + l)j Asu(2)LUmn, k.a)Asu(2)RUmn, kn) , 
k-i...k^ i—1 



^For the relation with General Relativity see 1661 1671 : for a recent discussion on the appropriateness of taking 
simple representations on nodes see | 18l 1191 l38l 139 40 41' and the discussion section at the end of this paper. 
^It can be compared to the analogous sign factor appearing in [7] and in I68| for gravity in three dimensions. 
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It can be written in diagrammatic notation as foUowfH 

1. jir, ,.. jiT, 



n ( )*o I 

ki...k5 4=1 J;\ 



k-i k2 

The occurrence of Imn/'^ as argument of ^bc in CH) can be understood as a consequence of the 
fact that the generator of SO{3) spatial rotations L can be written in terms of the generators Jl 
and Jr of SO{4) « SU{2)l(S) SU{'2.)r as L = + Jr. While at this level the specific form ^ is 
simply postulated, at the end of this paper we will give a check of it in a semiclassical regime. 

2.3 The boundary state 

Now we shift attention to the boundary state. As the Barrett-Crane model gives dynamics only 
to the integer spin subspace of Tips, we restrict to it. A spinnetwork state \T^,lran,in) € '^Fs 
describes a quantum geometry on a manifold . Such quantum geometry consists of 5 chunks of 
space (one for each node of the graph Fs), each meeting the other four chunks so that they identify 
in the whole 10 patches as prescribed by the connectivity the graph T^. This picture comes from 
the fact that the state jFa, Imm in) G is simultaneously an eigenstate of the volume operator of 
a region containing a node of Fs and an eigenstate of the area operator of a surface cut by a link 

of F5 [SSlEiElElj. 

Here, we are interested in a state iFs,;/) on Ti-r^ which is peaked both on the intrinsic and on 
the extrinsic geometry of . As such, we are looking for a (mildly) semiclassical state: it has the 
same connectivity of the state described above, but it is peaked both on area and on its conjugate 
momentum and similarly for volume. We take the following ansatz: 

1^5' 9) = X!X!^'^(*")*'o,0o(^mn)|r5,^mn,«n) (16) 

ill. 

with ^';o,0o(^nm) given by 

,T, fi \ ( ^ ST^ ST^ {Imn ~ lrnn){lpq — lpq)\ . ,(0) 7 /-, .7\ 

^"(0,00 ilmn) = exp - - 2^ 2^ a(„,„)(pg) , j exp (/-Ln^mn . (17) 

m<np<q y Imn Ipq m<n 

The matrix a(^rnn){pq) is assumed to be real and positive definite in order to suppress large fluctua- 
tions of Imn from the mean value Imli- As the quantity (j)mn appears only in the phase of (|17p and 
Ijnn is an integer, we can assume without loss of generality that it is in the range < (j>mn < 27r. 
It plays the role of mean value of the conjugate momentum to the representation Imn- To simplify 
the analysis, we restrict attention to a symmetric situation with Imli = lo and (pmli = 4>o- As stated 

closed loop gives the (super-)dimension of the representation space = (— 1)^^ (2j + 1), i.e. the dimension 

or minus the dimension depending on the representation being integer or half-integer. 




9 



before, we assume Iq ^ 1. Hence, due to the gaussian form of the state which is peaked on the 
value ^0 with dispersion -^/lo, we have that Imn is essentially restricted to be in the range 

(l-4=)^o<U. <(l + 4=)Zo. (18) 
v'o v'o 

This peakedness property is a kinematical property; if this state is to be considered semiclassical 
also in the dynamical sense, it depends on the specific spinfoam model for the dynamics chosen. 
In the following we discuss this issue for the Barrett-Crane model and show that the request of 
semiclassicality at a dynamical level fixes the angle (f>o to a specific value. 

It turns out to be useful to introduce the fluctuation Slmn and the relative fluctuation kmn 
defined as 

Imn — lo + Slmn = (1 + fcmn) ^0 (19) 

with 

Slmn 0( 1 (20) 



1 1 I'll' J — V fj — > 

The idea of taking coefficients with the simple form of a gaussian times a phase for the Imn part of 
the semiclassical state may appear naive, but it turns out to be an appropriate choice. We expect 
the form (|17p to be quite general, at least for large Iq, and not specific to the elementary graph 
Fs at hand here. Moreover it allows a parametrization of the corrections to the leading behaviour 
in terms of a polynomial multiplying the gaussian. The polynomial is a function of the relative 
fluctuation kmn with coefficients to be fixed. 

As the Barrett-Crane model is degenerate on intertwiner space, i.e. two states with the same 
representations Imn on links but different intertwiners at nodes evolve in the same way under the 
dynamics ^ defined by the model, we do not discuss the specific form of the coefficients /(««) 
and simply assume that /(*n) = 1- The coefficient C in has the role of normalization 
constant. The issue of peakedness on volume and appropriateness of the Barrett-Crane model will 
be discussed in section [H 

Due to the peakedness properties of the boundary semiclassical state, the dominant contribution 
from the Barrett-Crane vertex amplitude Wv(lmn,in) to the correlation function ^ comes from 
the range p8)) . In the following we give a detailed study of this contribution. We shall present it in 
terms of the function ^bc (jmn) defined in (jl5p in order to make the comparison with the existing 
literature easier. 



3 Integral formulae for the Barrett-Crane vertex 

In this section we give an integral formula for the Barrett-Crane vertex function defined in (|15p . 
We follow the analysis of [22] and reproduce the derivation here in order to fix the notation. Using 
the following formulgQ 

^In formula 1211 1 the diagrammatic notation f,^^^/ (h) = m — - — CX^) is used for the group element h in 

representation j. 
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> ■: -.k 



mi^j ji^m[ nil — CT>- 

- d^l{h) 



M^) ,3 ^ , (21) 

^ „ ^ ^ m.i — 



1714/ ?4 J-l\|7!4 '"4 -^—Ch^ "»4 

at each of the five nodes of equation (|15|) , we can express the Barrett-Crane vertex function as an 
integral over SU{2f [H] 



^BC(J12,..,J45) = (-1)^'^"" / ^^dy.[hu) n X^^-KKuhr^^) . (22) 

J5[/(2)= i<i,<5 l<m<n<5 



(there is a —1 factor for each half-integer loop; notice that this factor is cancelled by the analogous 
sign factor introduced in (|12p ). Now we rewrite it as an integral over (S^)^. A group element 
hk £ SU{2) can be described using the spherical coordinates [ipk,Qki4>k) on the group manifold of 
SU{2), i.e. on S^. Equivalently, hk can be described giving a unit-vector Vk in M"': 

Vk = ( sin -ipk sin 9k sin (f>k , sin ipk sin 6*4; cos 0^ , sin ■(/'fc cos 9k , cos V";!;) , (23) 

with the domain for {^pk, ^fc, 0fc) chosen in the following way 

0<1pk<Tr , 0<9k<'!T , 0<(t)k<2TT. (24) 

The Haar measure dfi{hk) on SU{2) can be given in terms of the Lebesgue measure on 

dfl{'il!k,0k,4>k) = sin^ iJkSm9kdipkd9kd(l)k ■ (25) 

As the S''^-solid angle is jdVl — 2tt'^ , the measure can be normalized to 1 introducing a factor 
(27r^)^^ in front of it. 

Using Weyl representation formula, the character in representation j of hmh^^ appearing in 
((22|) can be written in terms of the angle between the vectors Vm and Vn, 

{Vm, Vn) = COS $mn , (26) 

so that 

X^^-Hhrnhr^= ''''^''-^'^^- . (27) 
Sm ^rnn 

The integrand in equation (j22p depends only on scalar products {vm, Vn) of vectors in R^, hence it 
is S'0(4)-invariant. Using this invariance, we can choose coordinates on such that vi points to 
the "north pole" , V2 lies on the "Greenwich meridian" with latitude ■02 and 1)3 lies on the "2dim- 
Greenwich meridian" with 2-latitude ^3, ^3. In this way, we have fixed SO{A) invariance. With this 
choice of coordinates on S^, the five group elements ft-i, . . , /15 are described by the following vectors 

i;i =(0,0,0,1), (28) 

V2 =(0,0, sin?/'2, C0S?/'2), 

V3 =(0, sin -03 sin 03, sin -03 cos 613 , cos -03), 

V4 sin -04 sin 94 sin 04 , sin -04 sin 64 cos 04 , sin 1^4 cos ^4 , cos V'4) , 
115 = ( sin -05 sin 9 5 sin 05 , sin ip^ sin 9 5 cos 05 , sin ip^ cos 9 5 , cos "05 ) , 
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or equivalently by the nine angles u^, with i = I . . .9 

U = (-02, V'3,6'3,'04,6'4>4,V'5,6'5,<^5)- (29) 

Formula can be written as an integral over (S^)^ 

AC ■ \ ( ^^5:2J™„ f TT TT sin(2j„„» + l)$rrm .„„^ 

With the coordinates chosen above, the integrand does not depend on the angles "01; ^i? ^i; ^2, (/>2 
and 03, so we can integrate over them straightforwardly obtaining a IGtt^ factor: 

^Bc(ji2, ■ ■ ■ , j45) ^(-l)^^J"'" ^^^^^^ Jdilj2sm^ip2 Jdi/js y d6'3sin2 V3sin6'3 x 

X / dV'4 / dd^ / d(/)4 sin ■04 sin 64 / (i?/;5 / dd^ / d^s sin ^/is sin ^5 x 
Jo Jo Jo Jo Jo Jo 

^ TT Sin(2j„„i + l)$,n„(-02,V'3, 613, -04, 6*4, 04,05, 6'5,05) 
sin$„„(0J2,'03,6'3,'04,6'4,04,'05,6l5,05) 

The ten angles $mn between the vectors w„, v„i in can be written explicitly in terms of the nine 
angles uf. 

$12 = -02 , $13 = ^3 , $14 = -04 , $15 = ^5 , (32) 
$23 = cos^"^ ( COS 02 COS 03 + COS 613 sin 0^2 sin 1/13) , 
$24 = cos""'^ ( COS 02 COS 04 + COS 64 sin 02 sin -04) , 
$25 = cos^"^ ( COS 02 COS 05 + COS 65 sin 02 sin 05) , 

$34 = cos^"'^ ( COS 03 COS 04 + (cos 9^ 008 04 + COS 04 sin 6*3 sin04) sin 03 sin 04) , 

$35 — cos^^ ( COS 03 COS 05 + (cos ^3 COS 6*5 + COS 05 sin 6*3 sin 05 ) sin 03 sin 05) , 

$45 =cos^"^ ( COS 04 COS 05 + (cos ^4 COS 6*5 + cos(05 — 04) sin6'4 sin 6*5) sin 0^4 sin 05) . 

Calling D the 9-dimensional angular domain for the angles u^, formula (|3ip can be written in 
the following way 

Abc(ji2, . . . , J45) = (-l)^2j'„„ f J-J -Q ^-^ ^(2j;„„ + l)$„„(u)) (33) 

^ i=1...9 m<n 

where f{u) is given by 

, , IBtt** (sinwi)^ (sinM2)^ sinu3 (sinu4)^ sinu5 (sinwy)^ sinug . , 

^^"^^(2^ ni<™<„<5«in$™n(«) ■ ^ ' 

Moreover, the product J| sin ((2_;/m„ + l)$,„„(u)) in formula ([55]) can be written as a sum of 2^" 
terms in the following way 

n sin ((2j;„„ + 1)$™„(«)) = -1^ n (g+.(2,„„+i)*„„(..) _ ,-.(2,„„+i)*„„(«) 

l<m<n<5 7n<n 



2W E (-1)^"'"" exp (i J2 (-l)'-(2jmn + l)^mn{u)) , (35) 



6=0 
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where we have introduced the integer 6 = 0,.., 1023 and the binary digit notation 

b (612, ^13,^14, 615,^23,^24,625, ^34,^35,^45), (36) 

with b„in G {0,1} so that 

6 = ^ (612, 613, 614, 615, 623, 624, 625, 634, 635, 645) = (0, 0, 0, 0, 0, 0, 0, 0, 0, 0) 
6=1 ^ (612, 613, 614, 615, 623, 624, 625, 634, 635, 645) = (1, 0, 0, 0, 0, 0, 0, 0, 0, 0) 
6 = 2 ^ (612,613,614,615,623,624,625,634,635,645) = (0,1,0,0,0,0,0,0,0,0) 
6 = 3 (612,613,614,615,623,624,625,634,635,645) = (1, 1,0,0,0,0,0,0,0,0) 

6 = 1023 (612,613,614,615,623,624,625,634,635,645) = (1,1,1,1,1,1,1,1,1,1) 
Similarly, the (— l)2i;2j,„„ factor can be written as 

(_1)E2.„„ = cxp ( - ^E™<J-1)''""(2J™« + 1)^). 

Notice that, while the integrand in ([55)) is well defined for configurations u; such that <l>,„„(zi) = 
for some m < n, the function f{u) is singular at such points. These points correspond to degenerate 
configurations of the five vectors w„, with two or more of them coinciding. It turns out to be useful 
to introduce a quantity A-BCeUmn) defined as the integral on a domain such that a ball of radius 
e around degenerate configurations has been excised from D. The original formula for the vertex 
can be obtained taking the limit e ~* 0. 

Using formula ([55)1 and formula the Barrett-Crane vertex function can be written as (the 
hmit oQ) a sum of 2^° terms A'qI^: 

210-1 

Abc(J12, . . . , J45) = 1™ E (-1)^''""41(J12, • ■ -,^45) (37) 

b=0 

with 

4ce(ii2, • • -,^45) = / H du, f{u) exp {i E (-l)'""(2j™„ + l)($™„(ii) + tt)) (38) 

i=1...9 m<n 

with f{u) defined in ([M]) and $„i„(w) in (gH) , (021) ■ Formulae and for the Barrett-Crane 
vertex function will have a major role in the analysis of the following two sections. 

4 Dominant contribution from the Barrett- Crane vertex 

Thanks to the form of the boundary state q{lmm'in) discussed in section [2.31 we have that the 
correlation function ([5]) can be written in the following way 



W_dkmn Wv{{l + kmn)h) (1 + fcm'n') (1 + fc„i"„") ^';(,,0(,((1 + fcmn)^o) 
]Jdfc™„ H4((l + fcmn)^o) *io,0o((l + k„in)lo) 



^Notice that in II37I I. while the hmit of the sum is well defined, the sum of the limits is not. This amounts to take 
the Cauchy principal value of the addenda given by l|38|l . 
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where we have used the fact that, for ?o 3> 1, the sum over can be approximated by an integral 
over kmn and the eigenvalues of the area can be written as Amn ~ 87rGjv7^o(l + kmn + 0(fc^)). 

Using formula ([57]) for the Barrett-Crane vertex function, we have that the calculation of cor- 
relations of areas can be reduced to computing the following quantity 

Y[ dk^„P{kmn) S'^- 4c.((l + kmn)lo/2) *io,0o((l + kmnW (39) 

with P{kmn) given by a polynomial in kmn- 

In this section we use the integral formula (|38p for the Barrett-Crane vertex function to identify 
the dominant contribution to the quantity (|39p . Substituting the expression (|38p and the semi- 
classical state (|17p written in terms of the relative fluctuation kmn and of the mean spin la in the 
quantity (|39|) . we have 

I n dkmnP{kmn) (-1)'°^^"" A^^l^{{\ + fc™„)/o/2) -f kmn)lo) = (40) 

= / ndu./(.)e^5:(-)'-*'""H/"^(^^^^""*-^"^-^'') / n^fc™"^(fc™ox 

X C(Zo) exp ( - ^^0 X! '^^[m,n){pq)kmnkp^ exp Uq ^ (^(- 1)^'"" (w) - (/)o)fcmn 



m<n p<q 

9 



xexp(-izo E E"r™«)w((-i)'""'^'-"(")-'^-0(^^i)''''*^'(")~^^^^^ 

ni<n p<q 

where, in the first equality, we have exchanged the integral over kmn with the integral over the 
angles Ui while, in the second equality, we have performed the integral over kmn- The resulting 
expression has to be studied in the large Iq case. It turns out to be useful to analyse it in the limit 
^0 ^ oo the first: 

• The limit Iq oo and the regular configuration. In the limit we have that the gaussian 
appearing in the last line of (1^0]) becomes a product of delta functions 

l[s(^{-l)'^-"^mn{u)-cbo) - (41) 

As < ^mn{u) < IT and < 4>o < 27r, the delta function is non-zero only for (— l)**™" = -(-1 
Vm < n. Thus, in the limit — oo, only the term with 6 = in the sum p7p can contribute 
to the correlation function. Restricting to 5 0, we have to look for the angles Ui which solve 
the equation 

$mn(w) = 0o ym<n. (42) 
We find that there are only two solutions. We call them u° and Puq- The solution vP is given 

by 

U° = (l/'2,l/'3, 6*3, "04, 6*4, (?!'4,l/'5, 6*5,05) , (43) 
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with 



■02 =03 = 0^4 = 05 COS ^(-^) 
1, 



:27r- cos"i(--) 



cos-i(-l). 



This solution corresponds to the following configuration of the five vectors v. 

^ l^ _ 



VI =(0,0,0,1) 
/572 



W4 =(-- 



V2 = (0, 
/576 



4 ' 
1 



-I) 



^^3 = (0, V5/6, 



) : = (- 



/572 75/6 



J' 

5/3 



2' 2 4' 4' "'2 2 

These vectors have a 4-dimensional span and the angle between them is 



^mn = cos '^{Vm,Vn) = COS ^(-| ) 



(44) 



(45) 



(46) 



Hence we have that, in order for the boundary state to have non-zero correlation functions, 
the angle 0o in the semiclassical state ([T7)) has to be fixed to the value cos~i(— 1/4), at least 
in the limit Iq ^ 00. 

The other solution of (|^^ . P-o", corresponds to a second set of vectors Vn with a 4-dimcnsional 
span and can be obtained from the previous one by a reflection with respect to the 3-plane 
identified by the vectors vi, V2 and V3. It corresponds to exchanging with v^, or equiva- 
lently to exchanging 04 with 05. 



• Large lo and the dominant contribution. For large (but finite) Iq we have that, to the integral 
over u, only the u such that 

|(-l)^-$,„„(^i)-0o|<4= (47) 

V '0 

contribute to the area correlations. As the derivative of ^mn{u) at vP and at P^f' is of order 
one, we have that only the term 6 = can contribute and only the configuration u such that 
it belongs to a ball B^o centered in vP and of radius l/^/kj or to Bp^o (defined in a similar 
way). 

Taking into account the analysis presented above, we decompose the angular domain in the 
following way 

D, = Bao U Bp^o U 7^e . (48) 
Noticing that, as the integrand in psp is invariant under exchange of 04 with 05, we have 



and 



^ 2"'-l 
ABcijmn) = --^(2A'^lg^_^Jjrr,n) + lira ^ (- 1)^ ^BCe 0™" ) 

b=l 
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In terms of this decomposition, we have the following result: the terms ^bc^ (jmn) and A^^^ (j 
contribute to correlations only in a exponentially suppressed way: 

n dkr^nPiKnn) ^k'c°^((l + krnn)lo/2) *;o,0o((l + kmn)lo) = o{l^^) VA^ > , (50) 

m<n 

J(l + fc™„)?o/2)«'io,0o((l + fc„„O?o) -o(Zo^) ViV>0. (51) 

On the other hand, the term ^gpg ^ {jmn) contributes to correlations of areas giving a contribution 
which is suppressed only as a certain power of Zq: 

j n dkranP{kmn)A^l^^^{{l + fc„„0^o/2) *io,0o((l + fc™0/o) - 0(/o") (52) 

m<n 

for some fi to be determined. 

5 Large jo asymptotics of the dominant contribution: 
stationary phase approximation 

As discussed in the previous section, the dominant contribution to area correlation functions comes 
from the terms ^bcB-o^-^™") ^'^^ ^bcb _g{imn) of the Barrett-Crane vertex function, and the 
contributions they give are identical. Here we give a detailed analysis of such contribution studying 
the asymptoticf^ of ^g^g ^((1 + kmn)jo) for large jo, with fixed We recall that, thanks to 

the form of the boundary semiclassical state, we have 



mn ) 



JO 



< (53) 



In the following we view the relative fluctuation kmn as a small parameter, independent from Jq. 

The large jo asymptotic behavior of ^Bce ^ ((1 + kmn)ja) can be found using the method of the 
stationary phase [69]: the large A behaviour of a function F{X, k) of the form 



with f{u) smooth and with compact support in D and with Qk{u) smooth and having a single 
isolated stationary point u in the interior of D, is given by the asymptotic expansion 



^^(A, = (^)' e ^^^?'^(^) ( ^ «„(.) A-« + o(A-^)) . (54) 

The coefficients a„ (u) of the asymptotic series are independent from A and are determined in terms 
of derivatives of j(u) and Qk,{u) at the stationary point u. As a consequence, they depend on the 



'For a detailed discussion of the asymptotics of ^b^tj and of J^-qq^ see |22j and |23 ) . |24 | . 



16 



parameter k. In particular the zero order coefficient oq is given by 



f{u) e^^i 

ao(w) = —r==== (55) 



/|det(#S^r)l 



where the sign in ±ij is given by the signature of the Hessian g„ . g„ . | ■ For a Feynman-diagram 
algorithm for the coefficients a„(u) see for instance section 5 of [7D]. If the function /(u) does not 
vanish on the boundary of the domain Z), a next-to-leading-order contribution appears. It can be 
determined by integration by parts: if the phase Qniudo) has a stationary point within dD, an 
extra contribution B{X, k) — O(A^^A^^t^) = O(A^t^) appears on the right hand side of equation 

An asymptotic analysis of ^Bce- Unm) is possible thanks to formula 

(56) 



i=1...9 



5.1 The perturbative action 

To obtain the large jo asymptc 
of the phase Qk{u) defined as 



To obtain the large jo asymptotics of ^Bcg_ ((1 + kmn)jo), we have to study the stationary points 



Qk{u) - XI (1 + k^n){^mn{u) + Tt) , (57) 



m<n 

i.e., fixed the fluctuations /c^n^ we have to look for angles u.i G B^o such that 

dQk 



= 



dui 



(58) 



The function Qk{u) will play a key role in the determination of a perturbative action Sig{6lmn) as 
we shall show in next section. 

As noticed in [50], equation ([55]) has the same form of the Schlafli differential identity for a 
single 4-simplex having faces of area (1 -I- fc,„„)ylo (for some given reference scale Aq) and dihedral 
angles ^mn- Here we do not make use of this fact and look for an explicit expressions for <I',„„(u) 
and Qk{u) as functions of the parameter kmn- 

Our strategy is the following: we start studying the problem for kmn — 0, i.e. we look for angles 



uf^ in Bun such that 



dui 



(59) 



One can check that w° determined in the previous section as a solution of (|42)) is an isolated 
stationary point, i.e. it is a solution of equation (|59p. Moreover it is the only stationary point 
in Sgo, i.e. in a ball of radius I/a/jo = 0(fc) around vP . Then, we study perturbation theory in 
kmn ^ 1 around the stationary point. 



17 



Given a solution of (|59p . the stationary configuration u of equation (|58p can be determined 
perturbatively in fcmn- We write a solution of equation ([58| as a series 



S = + Sd) + 1^(2) + i^l(3) + . . . (60) 

with u^") = 0(fc"). In particular the linear perturbation is first order in fc,„„. Adopting the 
notation ™" (u) = <i>m„,i(u), we have that u has to satisfy 

= ^ (1 + fc™„) $™„,.(?2° + + is(') + . . .) (61) 

m<n 

order by order. Hence, expanding up to 0{k^) and recalling that — 0(1), u*-^-* — 0{k), m*^^-' = 
0(fc2), we find that vP, and m^^-' have to satisfy 

= ^ $m„,z(M°), (62) 
= kr,^n^mnA^°) + E ( E (^°)) ^ : (63) 

= J E ( E *™",.,i(u°))4')u|'' + E ( E ^mn$mn,^, (u°))^if + 



2 

j I m<n j m<.n 



iE(E*™-.(^"))"r- (64) 



2 

j m<n 

We are interested in the value of the phase Qk{u) evaluated at a stationary point u. Expanding up 
to third order we find 

g,(z2" + sW + ifi(2)^^-(3)_^ ^ ^ )^ 

= E (*m„(M°) +7r) + ( J2 fcmn($mn(S°) +7r) + ^ ( ^ $™„,^ (S°)) ^ ) + 

m<n 7n<n i m<n 

+ ^ f E ( E <^^n,,{u'))u'.'^uf^ + 2E ( E ^™n'i>".n.("°))"r^+ 
+E ( E $™n,.(s°))"f') + E ( E <i>™«...(^°))^!^^^f^^[^^+ 

i m<n / \ ijl m<n 



^E ( E ^™n$„m.»,(M°))zif^flf^ +3^ ( ^ <^^nM^°) 

ij 7n<n i j m<.n 

3E ( E ^mn1>™„.(S°))?2f) + ^ ( ^mnA^°)yfA+0{k^) 



Substituting in the above perturbative expression for Qk{u) the conditions (|62j|64p on 
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we have 



7(1) 



1 (2) 1 



7(3) 



2" '3! 



u^-' + ...) = 



(65) 



m<n 



ij m<n 



Notice that the dependence on m(^) appears only startmg from the forth order. This is a straight- 
forward consequence of the function Qk{u) being evaluated at a stationary point. Hence, to obtain 
the value of Qk[u) up to 0{k^) we need to compute only the correction vS^^ to the stationary 
configuration 

The correction u(i) can be determined in terms of vP and of kmn using the linear equation (|63p . 
Defining the matrix My = X]p<q '^pq,ij{^^) ^^id the "vector" N(^^n)i = ^7nn,i{u°), we have that the 
previous equation ([63|) is solved by 



{mn)i kmn with M^(mn)j 



(66) 



This gives the stationary configuration u to first order in fc„ 



/ cos-l(-|) \ 
cos-i(-|) 
cos-i(-|) 
cos-i(-l) 
cos-i(-i) 
27r-cos-i(-|) 
cos-i(-|) 
cos-i(-|) 

V cos-i(-l) y 



/ V^(_18fci2 + 7(fcl3 + fcl4 + fcl5 + fc23 + fc24 + fos) - 8(A:34 - ^35 - k4s)) 
(-18fci3 + 7(fci2 + fcl4 + fcl5 + ^23 + ^34 + ^35) - 8(A:24 - ^25 - ^45)) 



8 

/3/5 
8 

-^{kl2 + kis ~ ki4 ~ fcis + 2(fc24 + k25 + ^34 + ^35) " 4(fe23 + ^45)) 



g— {-18fcl4 + 7(fci2 + fcl3 + ki5 + fc24 + ^34 + ^45) - 8(A:23 - ^25 - ^35)) 
:^{ki2 - fcl3 + ki4 - fei5 + 2(fc23 + k25 + A:34 + ^45) - 4(fc24 + ^35)) 



(A;i3 + ki4 + k23 + k24 + ^35 + A;45 - 2{fci5 + fc25 + ^34)) 



g— {-18fcl5 + 7(fcl2 + fel3 + fel4 + fc25 + ^35 + ^45) - 8(A;23 + ^24 + ^34)) 
l^(kl2 - fcl3 - ki4 + fci5 + 2(fc23 + ^24 + ^35 + ^45) - 4(fc25 + ^34)) 



V 



(A;i3 + fei5 + A;23 + fes + ^34 + k45 - 2(fei4 + A:24 + ^35)) 



-0(fe2). 



The stationary configuration to first order in kmn is all what we need to compute a Taylor 
expansion in kmn of the phase Qkiu{kmn)) up to third order. Substituting (|66p in (|65p . we find 



Qki^'^ikmn)) — ^" ^ ^ ^mnkmn ^" ^ ^ ^ -^{mn)(pq)kmnkpq~\~ 



(67) 



3! 



m<n p<g r<s 
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with 





= Y1 (*mn("°)+^) -10(cOS-l(-l/4) 

= «'m«(w°) +7r = cos-i(-l/4) + 7r, 


+ 7r) , 


(68) 






(69) 


-^{rnn) {pq) 


ij u<Cv 




(70) 


^{mn) (pq) (rs) 


+ 3^ $™„,y(u°)W^(p,),W^(^,)j- . 


(71) 



i j I u<v i j 

The coefficients Qq and Bmn are easy to determine due to the fact that 

•/■O = *™n("") = C0S-l(-l/4) , 



(72) 



as shown in equation (pS]). The coefficients i^(mn)(p(j) and I{mn){pq){rs) can be calculated exphcitly. 



For Ki^rnn){pq) wc havc 



(mn)(p<j) 



/-9 


7/2 


7/2 


7/2 


7/2 


7/2 


7/2 


-4 


-4 


-4 


7/2 


-9 


7/2 


7/2 


7/2 


-4 


-4 


7/2 


7/2 


-4 


7/2 


7/2 


-9 


7/2 


-4 


7/2 


-4 


7/2 


-4 


7/2 


7/2 


7/2 


7/2 


-9 


-4 


-4 


7/2 


-4 


7/2 


7/2 


7/2 


7/2 


-4 


-4 


-9 


7/2 


7/2 


7/2 


7/2 


-4 


7/2 


-4 


7/2 


-4 


7/2 


-9 


7/2 


7/2 


-4 


7/2 


7/2 


-4 


-4 


7/2 


7/2 


7/2 


-9 


-4 


7/2 


7/2 


-4 


7/2 


7/2 


-4 


7/2 


7/2 


-4 


-9 


7/2 


7/2 


-4 


7/2 


-4 


7/2 


7/2 


-4 


7/2 


7/2 


-9 


7/2 


U4 


-4 


7/2 


7/2 


-4 


7/2 


7/2 


7/2 


7/2 


-9 



(73) 



where the following ordering for the fluctuations has been used, 

{kmn} = (^12, fcl3, hi, ^15, fes, ^24, ^25, ^34, ^35, ^45) • 

As expected from the symmetries of the boundary state and the fact that two links belonging to 
the graph can be either coincident, or touching at a node or being disjoint, we have that giving 
the following three components of the matrix if(mn)(pg) 



(12)(12) 



(74) 



9/3 7/3 /3 

"4V5 ' ^(i2)(i3) = +gy5 , i^(i2)(34) = 5 . 

is enough to reconstruct its structure. The matrix K(^„in){pq) has five negative degenerate eigen- 
values A_ = — \/l5, four positive degenerate eigenvalues A+ = \/l5/8, and a zero eigenvalue Aq 
corresponding to the eigenvector vo = 1/vTO (1, 1,1,1,1,1,1,1,1,1). The presence of a zero eigen- 
value was to be expected due to the fact that a symmetric fluctuation kmn = gives simply a 
small shift of the large jo . 

Similarly, for the coefficients (pg)(rs) there are seven classes of elements: 



'(12)(12)(12) 



'(12)(23)(13) 



'(12)(13)(14) 



189 

"so" 

141 
"20" 
453 
160 



'(12)(12)(13) 



'(12)(23)(34) 



'(12)(23)(45) 



347 /3 

I60V 5 
39 Is 
'20V5 ' 

loV 5 ■ 



'(12)(12)(34) 



14 



(75) 
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All the other components of the tensor I[mn)(pq)(rs) can be obtained by symmetry arguments. 
5.2 The perturbative measure 

Now we shift attention to the term ao(u(fc)) in the asymptotic expansion (I54p . Its general expression 
is given in equation (|55p . While the analysis above deals with the perturbative action, the following 
is a derivation of the perturbative measure. 

We evaluate the determinant of the Hessian of Qk{u) at the stationary point u, perturbatively 
in kjnn- Let's call Hij and H'j'j the 9x9 matriceJ^ 

^ ^ 9^Qk o-(o) _ d^Qo 



u"+uy'->+... 



We have that H,j = H^f + H[f + 0(fc2), with i/.^^ = 0{k). The matrices and H[f can 

be computed explicitly using the results of the previous section and are given by the following 
expressions 

The matrix Hf^ has four negative non-degenerate eigenvalues and five positive non-degenerate ones. 

They are all of order one. Then we have that the signatures of Hij and of Hj^j'' coincide, as their 
difference is of order 0{k). The determinant of Hij can be computed to first order in k using the 
following formula: 

det H = det {Hq (1 + H^^Hi + 0{k^)) ) 

= (deti7o)(l + TiiH^^Hi)) + 0{k^) . 

An explicit calculation gives 



det 



duiduj 



512 /5 3 17 

'10 - -fcia-f— (fci3-f /cM-f fci5)-h (76) 



177147 V 3 ' 2 2 

+ ll(fc23 + k2A + k25 + k34 + fcss + kir,)) + 0{k^) . 

Similarly, the function f{u) defined in p4p can be evaluated perturbatively at the stationary point 
u and we find that it is given by 

- {k23 + k24 + k25 + k34 + ^35 + ^45)) ) + 0{k^) . 

The ratio appearing in equation (|55p is given by 



^''Do not mistake Hf^ for The first is the Hessian of Qoiu) with respect to the nine angle Ui at the 

stationary configuration ft". On the other hand, J''(mn)(p9) the Hessian of Qk{u) with respect to the ten variables 

^mn^ at krnn — 
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Remarkably, to first order, it is symmetric under exchange of the ten kmn- This fact was to be 
expected on general grounds and provides a consistency check of the calculation. 

The last term in (j56p which needs to be computed is exp i ^ $mn (u) . Using the techniques and 
the notation of the previous section we have 

$,„„(^2" + + 1^(2) + 1^1(3) + ...)= E <S>mniu°) + ^ E ( E + 

ijl m<n ij in<n 

Comparing with the definitions ([70)1 . ([TT]) and (|72p . one finds that 

mn) ipq) ^^mn ^pq~\~ ('^^) 

m<n m<n mKn p<q 

+ "^i" y ^ ^ ^ ^ ^X ^'^^{mn)(pq)(rs))^m7ikpqkrs + 



Having computed all the needed ingredients, we conclude this section giving the asymptotics of 
the quantity ^BCg_ ((1 + Knn)jo)'- we have that 

A'bU, ((1 + ^™„)jo) = (^) ' 7^ (1 - ^ E + O(fc^) + O{jo')) X (80) 

The functions Qk{u{k)) and ^ ^mn{u{k)) contributing to the phase are given in equation (|67p and 
— 5 

([75]). The leading contribution to the asymptotic expansion is O(jo The correction of order 
0{k) has been computed explicitly. Corrections of order 0{k^) and 0{jQ^) can be computed using 
the techniques introduced above: they come from the next-to-leading-order expansion of ao('S(fc)) 
and from the leading order of ai(u(fc)). The extra term -B(joi K^n) appearing in (|80p is due to the 
fact that the function f{u) given in equation p4p does not vanish on dBu„. Such term is of order 
0{jQ^) and contributes to semiclassical correlation functions only in an exponentially suppressed 
way as discussed in section |4l 



6 Two-and three-area correlations on a semiclassical state 

Now that the dominant contribution from the Barrett-Crane vertex function to correlation formula 
([5]) has been identified and computed explicitly, we are ready to analyse correlations on a semi- 
classical state. The key result of sections |4] and [5] can be stated as follows: the spinfoam vertex 
amplitude defined in equation ()14p has the following asymptotics for large lo and Slmn of order 

W,{lo + 5l,nn) = ( n (2(^0 + SU^) + 1)^/) (-1)5:^'™.Abc((1 + ^-j^)lo/2) 

= Affl,,{SUn)e'^'oiSlmn) + ji^i^ + si„,n) ■ (81) 
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The perturbative action Sig{Slrnn) is given bj0 



— ^ — — ^0 

■m<n m<n p<q 

m<.n p<q r<s ^ ^ 

while the perturbative measure ^i„{Slmn) is given by 

9 , ^ 5L 



= C^"--^ (i + (^/ - 1) E ^ + • • • ) (83) 



with the dots standing for terms of order 0((^^^)^) and 0{Iq^) which can be determined using 
the techniques of the previous section. 

The function R{lo + 6lmn) in equation ()81|) contributes only in a exponentially suppressed way 
to semiclassical correlation functions as discussed in section [4] (see in particular equations (I50|) and 

m)- 

As large area eigenvalues can be written as 

Amn = SnGNlVilo + 5lr,M + 5Unn + 1) = 8^Gjv7 ?o f 1 + + + . . . ) (84) 

it turns out to be useful to compute correlations of ilmnl^a on the state ^'r5,(j('mn) before computing 
area correlations ( Apq )q and ( A^^ Apg Aj,s )q- We have that the leading contribution to the 
one-spin and two-spin correlation functions are given by 

(%^), = + 0(lAo), (85) 

(^V^^ = ^(*^--)«.)+W^o). (86) 

To compute the three-spin correlation function at leading order, an improvement of the boundary 
semiclassical state is required. We introduce a perturbative parametrization of the correction needed 
in terms of a polynomial in Slmn/lo multiplying the 'free' semiclassical state: 

*V0o(/o + 5/™0 = (1 + ci E ^ + • ■ ■) exp ( - i E E ^^^^^Sl^nSlpq) X 

— to V z ^ — ^ ^ — ^ to ^ 

'm<n m<n p<q 

X exp -i(j>o E (^0 + Slmn) ■ (87) 



-'^-'^We have called perturbative action the function Sig{Slmn) due to its similarity to a (lattice) field theoretical 
action 1711 . Moreover in the following we use some standard techniques from perturbative quantum field theory, 
such as splitting the action in a free part and an interaction that can be treated perturbatively in the path integral. 
However the name action is not completely appropriate and can generate some confusion: an appropriate name for 
the function Sig{Slmn) is perturbative Hamilton function as it plays the same role the action evaluated on a classical 
solution plays in WKB semiclassical perturbation theory. 
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The leading order correction is parametrized by the single constant ci . The three-spin correlation 
function can be computed perturbatively viewing the cubic term in the second line of (j82p as an 
interaction term which perturbs the free quadratic action. Using Wick's theorem we find: 



"-O ^0 ''0 ''0 / 

By we mean the sum over the primed indices appearing in the formula, X]m'<n' instance. On 
the second line of (|88p nine terms appear corresponding to permutations of the indices of the term 
written explicitly. Similarly, on the third line three terms appear corresponding to permutations of 
the indices of the written term. Taking into account the fact that the order O{l/lo) correction to 
(123 is 

+ (Cl +Nf-^)J2 - + Oil/ll) . (89) 

we have that the joint cumulant for three spin fluctuations (i.e. the connected three-point correlation 
function) is given by the following simpler expression 

SlfYin ^^pq Sirs \ ( i ^^mn ^^pq \ / ^^rs \ , \ , n / ^^rtin \ i ^^pq \ i Sirs \ /r\r\\ 

{—, T^-r^i - [(—, T^)ii^)i + Pe™. j + 2(-^)5(-p)g(— )q = (90) 

<o to '0 ^ to '0 to ^ to to to 



Using the previous results and the expansion (j84p we have that the area expectation value and the 
two-area correlation function to leading order are given by 

{A„,n)q = 87rGAr7 lo (1 + Oil/lo)) , (91) 

/ (iK — a)7^ ., , s 

{Aran Ap,) , - {Amn)q{Apq)q = [SnGNlf ( -i™^ + 0{l/ll)) . (92) 

The three-area correlation function ( A„„i Apg Aj-s )q can be computed expanding up to order 0(1/ Iq) 
and 0(^/^„/^q) the expression ([51)) for Amn 

=87rGA.7;o(l + ^-^ + 0{l/ll) + (1 + 4 + 0{l/ll))^-^ + 

\ Z60 OtQ OtQ to 

1 /^^mn\2 „/l/5U 



+ (-^ + 0(l//o^))(^)' + 0(^(f)^)) 

OtQ to tg to ^ 
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' [r' s' )(rs) 



and using expressions (ISS)) and (|89p . We have that the three-area joint cumulant is simply given 
by the following expression 

i^mn ^pq -^rs)q ~ ( (^mn ^p(})q (^rs) g + Perm.) + 2 g (Apg) ^ (Ars) g = (93) 

= (8^Gjv7)^ (i ihm'n'){p'q'){r's'){iK - a)(;,i,„,)(„„) {iK - {iK 

7 Correlations in perturbative quantum Regge- calculus 

Our aim in this section is to compute correlations of geometric quantities within the framework of 
perturbative quantum Regge calculus. The motivation behind it is to compare such correlations 
with the ones computed from Loop Quantum Gravity in the preceding section. For an introduction 
to Regge calculus [1] and to the many issues in quantizing it, the reader can refer to the following 
recent reviews [7H [731 UM ■ Here we give only a brief introduction. The attitude we shall adopt is 
that to consider Regge gravity as ordinary General Relativity but with the metric g^^^ (x) restricted 
to belong to the class of piecewise-flat metrics |75[ , instead of belonging for instance to C°° . 

Let's consider a four dimensional manifold M with a Riemannian piecewise-flat metric gfj,u{x) 
|77| . A structure of this kind is called a simplicial manifold and is completely described by the 
connectivity of a finite number of vertices and by the l eng ths of edges connecting two of the vertices. 
Some nomenclature: we call C the simplicial complex^ describing the connectivity of the vertices, 
Co = {"^1, ■ ■ ■ ,vn} the ensemble of its vertices, Ci the ensemble of its edges, C2 the ensemble of its 
triangles, C3 the ensemble of its tetrahedra, C4 the ensemble of its 4-simplices. Moreover we call: 
Cij the edge from the vertex Vi to the vertex Vj and its length, tijh the triangle with vertices 
Vi,Vj,Vh and Aijh its area, Tijhk a tetrahedron and Vijhk its volume, T^ijuu a 4-simplex and Sijhki 
its 4-volume. Moreover, we call 9(ijh)ki the dihedral angle between the tetrahedra Tijhk and Tijhi 
which share the triangle tijh- 

As well-known, the Einstein-Hilbert action for a manifold with boundary is given by [781 179] 

S[9M = T^^ I d'^V9R+^^ [ d'xVhK. (94) 

For a piecewise-flat metric gjw'^"^"\x), we have that S[gj£'^"^"\x)] — S'Roggc [C, Lm„] [75l [76]. 
with the Regge action given by [?[ lED] 

'5'Roggo[C', Z/„i„] = - — — — ^ ] Aijh{Lmn)(^'2'^ ~ ^ ] 9(^ijh)kliLmn)^ ^ 

i<j<h k<l 

+ Aijh {Lmn) f ^ ~ 9{ijh)kl {Lmn) 

^ i<j<h k<l 

tijh&dC2 Tijhk ,Tij hi SidCs 

The quantum theory can be introduced a la Feynman '81] as a sum over piecewise-flat geometries. 
This involves a sum over both the connectivity and the lengths, with a specific measure and with 



-'^ There is a restriction a simplicial complex C has to satisfy in order to describe the connectivity of a simplicial 
manifold. The reader can refer to the lectures 1741 for details. 
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an appropriate gauge-fixing procedure in order to avoid multiple counting of a given piecewise-flat 
geometry. An attitude one can provisionally adopt is to freeze the connectivity to a given simplicial 
complex C , and integrate only over the lengths of the edges of C. Within this setting, the problem 
of the choice of the measure fici^mn) has been widely discussed in the literature. Following the 
point of view we have adhered here, such measure can in principle be determined starting from 
DeWitt measure for g/ji^' (a;) and computing the Faddeev-Popov determinant for the change of 
variables to the invariants L„,„ l82l [83l [84l [85l [86l : 



/ 



^igj^'^^-Hx)] SixigjS^'^-Hx))] Al^[gj£'^-Hx)] = J {U dL^n) • (95) 



The measure n^{Lmn) turns out to be highly non-trivial to compute^ and is expected to be given 
by a product over 4-simplices of functions depending only on the single 4-simplex edge-lengths. On 
the other hand, taking the point of view that the measure comes from a discretization of a formal 
measure in the continuum, a class of ultra-local measures fJ.^{Lmn) have been proposed and studied 
numericalljl^: they are given by 

I ( n ^■^"») ^^ciLran) = / ( H ^mn^imn) ©(^nm) (96) 

with the exponent N generally taken to be -1-1 or —1. The step function Q{Lmn) in the right 
hand side of (ITOj) ensures that the edge lengths satisfy the triangle inequalities and their higher 
dimensional analogs (see for instance [23] or appendix El. 

As always in a functional integral, boundary conditions have to be specified. Here we assume 
asymptotic flatness. Wc follow Rocek and Williams [351 ISZ] and study the theory perturbatively 
around a classical solution of Einstein-Regge equations. The classical solution we consider is flat 
space, i.e. a configuration which has zero deficit angle at each triangl^H- This amounts to find a flat 
configuration (C, L^^j), and quantize perturbatively the fluctuations SLmn = Lmn — -^mn around 
this background. The fluctuations are assumed to vanish asymptotically. In this setting, the 
problem of the choice of functional measure is avoided, at least at the lowest order in perturbation 
theory, as it turns out to be given by 

^^c{LL. + 5L^n)^^^c{LL.) (l+ E E E^(1)(P,)^§^---) (97) 

m<n ™" m<np<q ™" 

with specific coefficient^^ A*(mini) - (mfcnfc)- Perturbatively, one can study expectation values of 
products of localized operators, such as the correlation of volumes of two tetrahedra belonging to 



I'^For an analysis in two dimensions see [87] and 1881 1891 l90l 1911 . 
^^See for instance | 72l . |73l and references therein. See also |92| . 

^^In id the deficit angle at a triangle tijf^ is defined as Sijh = 2it — X]fc<i ^(ijh)kl- It can be shown that the 
densitized curvature for a piecewise-flat metric is related to the deficit angle in the following way; g{x)R(x) 

E ' 

^^Notice that we have assumed the zeroth order coefficient to be different from zero, /j''^) = AtcC^mn) ^- This 
input can be understood in the following way: it is needed in order to match with the continuum functional measure 
determined perturbatively. 



[2eijfi / cfia S^"^^ {x — x{a)) ), where x(a) is the embedding of the triangle i^jj, in A4. 
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{SVijhk SVi'j'h'k')Q — 



[ n dSL^n (l + E M^'^ !& + •••) 5V,,hk m^ww e*^c(^™" + '^^"») 

m<n 



J H dSLran (l + E /^^'^ Lt7 



This integral can be performed as follows. Choose a region which is a subset of with the 
topology of a 4-ball and such that the tetrahedra Tijhk and Tiiji^'k' belong to its boundary. Call 
5*3 its boundary. The action can be written as a sum over the two regions, S(j = Sb + S(j_g. Then 
the integrals in the correlation formula can be split in the following way: 

J Y[ dSLmn W{SL,nn) SVijhk 5Ve j> h' k' *o('^imn) 

{5Vijhk 5Vi:j:h'k')o = — "'^"'"T^'^^ (98) 

where we have defined the functions W{6Lmn) and ^Q{6Lmn) 

W[5L„^,.)=( X{d5L„^,^ (1+ mS(„.„)%^ + ...) e^^-(^"» + ^L„„0, (99) 

m < 71 m < n 

ra<.n 7n<n _ 

em„eint(,C-B} e„„eint(C-S) 

which depend only on the length-fluctuation of the edges belonging to 5*3. The function W{6Lmn) 
plays the role of transition amplitude kernel. To lowest order in perturbation theory, it is given by 
the exponential of i times the Hamilton function, i.e. the perturbative Regge action evaluated on 
a classical solution corresponding to the boundary condition (5*3, L^„ + 6Lmn)- To higher order 
a non-trivial measure multiplying the exponential is expected to appear. On the other hand, the 
function ^oi^Lmn) describes the state of the system. In the present case it codes the asymptotic 
flatness boundary condition. The name it deserves is 'perturbative vacuum state'. It can be eval- 
uated perturbatively, and to lowest order it is expected to be given by a gaussian peaked on the 
mean value {SLmn) = times a phase which codes the mean value of the momentum conjugate to 
the variable 5Lmn- As such, it is peaked both on the fluctuation of the intrinsic and of the extrinsic 
curvature of 6*3. 



Here we take as region a single 4-simplex S 12345. In this case we have that the vector of 
edge-length fluctuations has 10 components {Lmn} — ((5Li2, . . , 15L45). Moreover we assume that 
the background edge-lengths are all equal, = Lq^ i.e. we have as background the regular 

4-simplex. As declared above, our aim is to compute correlations of areas, volumes, angles, etc. 
in this setting and compare them with the ones computed at the vertex amplitude level in Loop 
Quantum Gravity in section [6l Even if not strictly needed, a possibility is to change variables in 
the functional integral and choose variables closer to the Loop Quantum Gravity ones: instead of 
edge-length fluctuations, as variables we choose fluctuations of the areas SAhki of the 10 triangles 
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thki ■ In this setting, the transition amphtude kernel is easily written: it is simply given by 



WiSAi23, ■ ■ , -^^345) - ^^{Ao + 6Ahki) e'^^-ssoi^o + SAhh) (iqi) 
with 5Rcggc(Ao + SAhki) given bj{3 

SRegsMo + SAhkl) ^ J2 {Ao + SAhkl){^n-e^hkl){Lo + SL^n{SArst))) ■ (102) 

^ l<h<k<l<5 

Notice that in the case of a single 4-simplex the Regge action and its Hamilton function coincide. 
Given the symmetries of the background, the function fi{AQ + SAfiki) to lowest order in perturbation 
theory is given by 

liiAo + 6Ahki) = 1 + J- ^^hM + --- (103) 



On the other hand, the boundary state 5'o(^^/!jci) is much more difficult to compute from first 
principles due to its dependence on the connectivity of C-B and on the asymptotic condition on 
the outer boundary. However one can parametrize its form: as discussed above it is expected to be 
given by 

T / r < \ /i '^1 f A \ 11 Cl(rst)(uvz)[Ao) 

An ^ — ' ZSnLrM ^ — ' An 

" l<h<k<l<5 l<r<s<t<5 ' 

l<u<v<z<5 

^ exp-i— i^tjhiAo) SAijh ■ (104) 
sttG-w — 

l<i<j<h<5 

with a(rst){uvz){Ao) and ipijhiAo) dimensionless. The quantity ^pijhiAo) codes the mean value of 
the momentum conjugate to (87rGAr)~^<5Ay7,. In order to describe flat space, it has to be chosen to 
its classical value which is given by the derivative with respect to SAijh of the Hamilton function 

(Una): 

dAi-jh 



^,jh{Ao) = SttGn ^T^' iAo) - TT - %^(Lo) . (105) 



The function 4'ijhiAo) is exactly the angle between the normals to the tetrahedra Tijhk and Tijhi 
which share the triangle tijt and belong to a regular 4-simplex, i.e. it is the Regge version of the 
extrinsic curvature. It can be easily evaluated (see appendix E]) and its value is cos~^(— 1/4). 



7.1 The perturbative area-Regge-calculus action 

In this subsection we determine the perturbative expansion of the area-Regge-calculus action for 
a single 4-simplex (|102p . We start discussing the background: as shown in [S^ using Bang's 
theorem, a non-degenerate 4-simplex which has all its ten triangle-areas equal is necessarily regular, 
i.e. it has all its edge-lengths equal. Then we need to compute the derivatives of S'R,eggc(^o + 
SAhki) with respect to the area fluctuations. Despite the fact that an explicit formula for L,„„ = 

^^In the case of a single 4-simplex, for each triangle tj^^ there are only two tetrahedra sharing it and hence only 
one dihedral angle 6(ijh)kh which we have called f?(jj7,). 
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Lmn{Ai23, ■ ■ j^345) IS missing, nevertheless the derivative of Lmn with respect to Ahki evaluated 
at ^0 = -^Lq can be easily found using the inverse function theorem: let's introduce the matrices 

U(hkl){mn){Lo) = CtT^^\ ^o) and V(^,nn)(hkl){M) = a ("^o) = [U~^) (^rnn)(hkl) (106) 

O-Umn OAhkl 

so that we have SL„in = J2h<k<i '^(mn)(hki)^Ahki (see appendix!^ for details). The first, second and 
third derivatives of the action (|102p with respect to area-fluctuations and evaluated at the regular 
configuration are given by the following expressions 

^^(^o) = ^{n-e,ki{Lo)) = -i^cos-i(-l) , (f07) 
dAhki SttGat SttGn 4 

^^S'rckrc . , X 1 d9hki 



(^o) = - %)(-*)(io) , (108) 



dAhkidArst SttGn ^ dL 



P^fl (^O) - - E E ^("m)(rst)(io) V(pq)i^uvz){Lo) X 



dAhkldArstdAuvz ^'^^^ ,n<np<q 

a2 



-(^o)-E E ^(^o)%)(aM(^o) (^^o) . (io9) 

i<j a<b<c ■' 



In equation (|107p we have used the fact that, as a special case of Schlafli differential identity, we 
have 



dLr, 

l<h<k<l<5 



= 0. (IfO) 



Using the explicit expressions for the first and second derivative of the area and of the dihedral 
angle with respect to the edge lengths and the matrix V(mn){hki){Ao) given in appendix [XI we have 
evaluated expressions (|108p and (|f 09p . Thanks to the symmetries of the background and to the fact 
that two triangles in a 4-simplex can be cither coincident, or sharing an edge, or sharing a single 
vertex, we have that the matrix (|f OSp is completely fixed once the following three components have 
been computecf^ 

9/3 „ I 7 



'Rcggc / ■ N 4 V 5 C'jSRoggc , , ^ V ■-' fTTT\ 

-[Aa) - , — [AQ)-—T^—-r, (111) 



9^1239^123' ' S^GjvAo' 9^1239^124 ' ' SttGnAo 



9^1239^145 ' ' St^Gn Aq ■ 
Similarly, in the case of the third derivatives, we have that the following seven classes of components 



The matrix 11081 1 has a zero eigenvalue with area-fluctuation eigenvector {SAi^^i } corresponding to a homoge- 
neous expansion (or contraction) SAf^j-i = SAq. Notice that such homogeneous expansion can be absorbed in the 
background value, Aq — > Aq +5Ao, which is assumed to be large. For an interpretation of this conformal zero mode 
as a gauge degree of freedom see [94] . 
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completely fix all the components of p09p . Introducing the dimensionless coefficients 

K(hki)(rst) = SttGat Aq — — ^^^—{Aq) , (113) 

I{hkl) (rst) (uvz) = SttGjv Al ^^A^^^^^A ^"^0 ^ ' ^ 

we have that the perturbative expansion in fluctuations of area of the single 4-simplex action is 
given by 

5Reggc(Ao+M„fcO=^^(V'0 {Ao + SA,,kl) + ^ E ^^^^^J^SA,,kl5Arst + 

^ ^ h<k<l h<k<lr<s<t ° 



+ E E '^^"^g"""^ 6A,k^6Ar.,6A^.. + O(^) ) . (115) 

■ ft<fe<; r<s<t it<t)<z 

7.2 Correlations of area fluctuations: comparison with the Loop Quan- 
tum Gravity calculation 

Within the framework introduced above, computing correlations of area fluctuations in perturbative 
Regge-calculus is straightforward. However, before discussing them, we would like to stress the 
following facts: 

- the value of the dihedral angle tpijh{Ao) deflned in (|105p . p07p coincides with the phase (j>o 
defined in (|17p and determined by equations ((4T|) and (|46)) : 

- the three numerical coefficients which characterize the matrix iir(/ifei)(rst) defined in (jll3p in 
terms of the second derivatives with respect to the area of the Regge action determined (llOSp 
and evaluated in (jllip coincide with the three coefficients determined on the Loop Quantum 
Gravity side in equation ([74| : 

- similarly the seven coefficients (|114p . (|109p . (|112p match with the seven coefficients ([75]). 
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Taking into account these facts, we have that the perturbative area-Regge-action 5'Rcggc(^o+^^/ifci) 
defined in (|115p and the function Sig{5lmn) introduced in ((82)) to describe a perturbative regime of 
Loop Quantum Gravity with the dynamics implemented using the Barrett-Crane model match, at 
least up to third order, provided that we identify the spins and the areas in the following way: 

OTTO- AT OTTG-jv 

On the other hand, within the Loop Quantum Gravity framework, the relation between spins and 
eigenvalues of the area operator is well-known [55l |56l [57j. Using formula ([84]) we find that, at 
least to lowest order, the relation between spins and areas coming from Loop Quantum Gravity 
kinematics and the relation we find here at a dynamical semiclassical level (jllSp are in accorcF^. 

In order to have that the correlations of areas computed in perturbative area-Regge-calculus 
match with the ones computed in section [5] from Loop Quantum Gravity, we need much more than 
a matching of actions: a matching of functional measures is needed. This can be accomplished 
choosing a specific perturbative measure on the Regge calculus side so that it reproduces the 
perturbative measure ([83]) . To lowest order, this amounts to the choice /ii — Nf — 55 in (jlOSp . 



8 Discussion 



The result presented above establishes a correspondence between the non-perturbative dynamics 
of Loop Quantum Gravity on a semiclassical state and perturbative-Regge-calculus on its vacuum 
state. The many assumptions behind this result are discussed in section [2] on the Loop Quantum 
Gravity side and in section [7] on the Regge calculus side. A fact we would like to stress here is that 
such correspondence is a non-trivial result, in the sense that it cannot be established on general 
grounds without working out the details as done above. Some steps in the analysis that a priori 
can go wrong - but actually turn out to work - are the following: 

• the correspondence has been shown to hold for a specific model for the dynamics in Loop 
Quantum Gravity - the Barrett-Crane model - but there is no guarantee that it keeps on 
holding for other models for the dynamics. Obstructions can arise at two levels: 

(i) at the non-perturbative level, due to a different dynamics in the 4-valcnt spin-networks 
sector. For instance, using the Hamiltonian constraint introduced by Thiemann in [951 
[5B][S7] we have that the matrix elements (O turn out to be zero. Therefore, no matching 
with Regge calculus is to be expected; 

(ii) at the 'semiclassical' level, due to a suppression of spin-spin correlations. Before [551 
[26] . the Barrett-Crane model was suspected to suffer from this pathology. Within the 
framework discussed in the present paper, this scenario can be stated as follows: there is 
no superposition of spin-network states with graph T5 as in ([T6|) such that the two-spin 
correlation function is different from zero and of order O{l/lo) as in (|86|) . In the case of 
the Barrett-Crane model we have found that this dynamical requirement fixes the phase 
00 in the boundary semiclassical state (see equation ([iTjl V At a following stage, we have 
compared such phase with the dihedral angle appearing in (|105p and found the two to 

^^Notice that the matching between the kinematical relation II84I I and the dynamical semiclassical one 1 11161 requires 
the Immirzi parameter 7 to be equal to one. This is expected to be due to the specific form of the Barrett-Crane 
model for the dynamics and not to be a general fact. 
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match. On the other hand, assuming the correspondence with Regge calculus, one could 
start with a superposition of spin- network states such that, at the kinematical level, it 
describes a state peaked on the intrinsic and the extrinsic 3-geometry of the boundary 
of a regular 4-simplex as done in [551 US] • Then, the calculation of spin-spin correlations 
comes as a test of semiclassicality at the dynamical level. If these large scale correlations 
are present for the recently proposed new models for the dynamics [TSl [TH [351 IMl HOI HI] 
is still an open issue. The same remark is expected to hold also when, instead of using 
the spinfoam formalism, the dynamics is implemented using canonical methods as for 
the graph changing Master constraint [55] . 

• in our analysis we have chosen an highly symmetric semiclassical state supported by a graph 

which is dual to the boundary of a 4-simplex. Provided that obstructions (i) and (ii) are 
avoided so that the dynamics turns out to be effectively described by a perturbative action 
like the one in equation (|82[) . a question one could ask is if the 10 x 10 matrix -R'(mn)(pg) 
and the 10 x 10 x 10 tensor I(rnn){pq){rs) in (|82[) are completely fixed by symmetry so that 
their matching with the ones appearing on the Regge calculus side in equation can be 

predicted a priori. The answer is clearly no as one can respect symmetry and still have 3 
free parameters at the level of if(mn)(p(3) and 7 free parameters at the level of I{mn){pq){rs)- 
The reason of the matching comes from the correspondence of the stationary phase equation 
([58)1 with the Schlafli differential identity for a single 4-simplex (|123p . Hence we expect on 
general grounds that the matching keeps on holding at the following orders in the Taylor 
expansion. What is non-trivial here is the matching of dimensionful parameters, namely the 
(BttGtv)"^ in front of the perturbative action. In fact Schlafli identity cannot fix the overall 
scale. The correct dimensions come from the fact that in ([5^ the dominant terms are of the 
form (5^m„//o^^, together with the relation ([6]) between spin and area eigenvalues confirmed 
in (|116p . On the other hand, the subleading contribution coming from the term ([79]) gives a 
correction which amounts to a shift — > Zq + 1- 

• in our analysis we have computed perturbatively the measure ^ii„{5lmn) and found that to 
lowest order it is given by a non- vanishing constant plus a linear term in Slmn/lo with a specific 
coefficient (see equation ([53])). While on the Regge calculus side there is some uncertainty 
on the choice of the appropriate measure, the fact that perturbatively it starts with a non- 
vanishing constant is certainly a requirement if one wishes to make contact with perturbative 
quantum gravity in the continuum (see for instance 92\ for a discussion and an example of 
pathological measure in quantum Regge calculus). 

Moreover, using the techniques introduced in sections[4]and[5l our analysis can be extended to higher 
orders in perturbation theory. A matter of interest here is the study of higher order corrections to 
the measure fJ.ig{Slmn)- A detailed analysis could clarify the following issues: 

- going at least to quadratic order in Slmn would allow to check if the perturbative measure 
coming from Loop Quantum Gravity matches with the expansion of the diffeomorphism- 
invariant measure (|95p /if_^;„j(^o + SAhki) determined on the quantum- Regge- calculus side 
following [82J and t83»i84, 85, 86], or not. We argue that something new about diffeomorphism 
invariance in Loop Quantum Gravity can be learned from this comparison; 

- the perturbative expansion of the coefficients a„(u(A;m„))Z(^" appearing in the expression (I54p 
for the asymptotic expansion can be exponentiated giving an 'effective action' Sfl^{Slrnn)- 
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On the other hand, the result (|116p provides a dictionary between spins and dimensionful 
quantities. By dimensional analysis we have that from ao can come only terms independent 
of Gjv, from ai only terms which scale as SttGn and similarly for higher orders. This suggests 
that the contribution to the perturbative action coming from the measure ^ii„{5lmn) can be 
written in terms of an higher-curvature Regge action [5511100] . with specific coefficients for the 
higher-curvature terms fixed by the non-perturbative theory. While suggestive, this conjecture 
needs a detailed analvsij^ of the expansion in fc,„„ of the an(u{kmn)) to be confirmed or 
rejected. 

These aspects clearly deserve further study. On the other hand, the analysis we have presented 
points out a difficulty in the matching of area-Regge-calculus variables with the genuine Loop Quan- 
tum Gravity degrees of freedom. The point is the following: the framework of perturbative-area- 
Regge calculus of section [7] allows to compute correlations of fluctuations of every geometrical object 
belonging to the 4-simplex S12345, and not only correlations of area-fluctuations. For instance, cor- 
relations of fluctuations of tetrahedron- volumes or of angle-between-triangles are straightforward 
to compute as such objects can be written in terms of fluctuations of areas (see for instance ap- 
pendix The same could be done on the Loop Quantum Gravity side, introducing appropriate 
functions of the spins in order to represent volumes and angles. Clearly, the correlations computed 
in this non-standard way will match with the ones computed in perturbative Regge-calculus by 
construction. The problem comes when one realizes that in Loop Quantum Gravity the volume 
operator [551 IH] [ST] and the angle operator jlOlj can be introduced at the kinematical level and 
their eigenvalues depend on the intertwiners at the nodes, and not simply on the spins on the links. 
As the Barrett-Crane model gives trivial dynamics to intertwiner space, we have that correlations 
of intertwiners will be unrelated to correlations of volumes and angles computed in perturbative 
quantum Regge-calculus, a fact already noticed in [26] and discussed in great detail in [102| . A 
possible reaction to this fact is to blame the Barrett-Crane model and to abandon it. Recently, 
some new models for the dynamics have been proposed [THl [THl [Ml [3H1 HOI HI] j all of them sharing 
the property of giving non-trivial dynamics to intertwiner space. From this point of view, a problem 
which urges to be addressed is to repeat our analysis for the new models, using a semiclassical state 
^v^^'^i'^mmin) given by a superposition both on spins and on intertwiners with the property of 
being peaked exactly on the intrinsic and the extrinsic 3-geometry of the boundary of a regular flat 
4-simplex, instead of the naive state 5'io,0o (^mn) of equation pT|) . Luckily, a state of this kind can 
be built thanks to the result [13], (see for instance [102] '). 

Computing correlations for the new models and comparing them with perturbative-Regge- 
calculus provides a non-trivial test of their viability as models for quantum gravity. Moreover 
it suggests a new way of interpreting the Barrett-Crane model W^c{lmn)'- namely to understand it 
as a 'partially integrated' model. Suppose for instance that a new model Wncw('mn, in) exists with 
the following property: the sum over intertwiners on the state ^'p~j!'™(^mn, *n) satisfies 

in 

When using such a model to compute only spin correlations the agreement with perturbative- 
Regge-calculus is guaranteed by our analysis of the Barrett-Crane model. On the other hand, when 
computing correlations of a function of the intertwiners such as the volume, the new model should 

^'^A preliminary analysis of the matching of the quadratic orders in Slmn of ag and ai with the expansions of the 
curvature-squared and the curvature-cube Regge action is under way. 
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satisfy the non-trivial requirement 

W^ncw('rrm, «n) Vi('lm, «l)V2(/2m, h) ^r's^'^C^mn, in) = W^Bc(^mn) Vl('mn)^2(^ mn) ^^0,00 (^^^^^) 

with the non-standard expression Vi{lmn) for the volume in terms of ten spins, an expression which 
can be checked against the area-Regge-calculus expression for the volume of a tetrahedron in terms 
of the ten areas of the 4-simplex it belongs to. While finding a model Wnow(^mn, in) which matches 
in this way with the Barrett-Crane model could be interesting, this is not at all a requirement for 
the new models [TSl [HI [351 [Ml [13 [S] • The weaker requirement they should satisfy is to have a 
partially integrated version with properties analogous to the ones postulated above for the Barrett- 
Crane model and such that it admits an effective description of the form (|8T|) . This amounts to 
changing the dynamics in the non-perturbative regime, while leaving the semiclassical description 
unchanged in form. 

There are three major open issues that the present work does not address: 

• the problem of going beyond the vertex amplitude level, namely the study of composition (a 
problem addressed at an exploratory level in [26]); 

• the problem of the sum over two-complexes (see [1] for a discussion), which corresponds to a 
sum over the connectivities on the Regge calculus side; 

• the relation with (continuum) perturbative quantum gravity on flat space, possibly viewed as 
an effective field theory description [1031 1104] . 

We hope that the analysis we have presented here could provide a safe starting point for the study 
of such problems. 
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A The geometry of a 4-simplex: length and area variables 

In this appendix we collect some formulas describing the geometry of a 4-simplex and determine 
the Jacobian of the transformation from variation of length- variables to variation of area- variables 
for the regular configuration. 

A 4-simplex is the convex hull of a set of 5 affinely independent points {ui, . . , us} in a Euclidean 
space of dimension 4. Subsets of {wi, . . ,^5} define 10 edges of length , 10 triangles tijh of 
area Aijh and 5 tetrahedra Tij^k of volume Vijhk- We call Sijhki the 4- volume of the 4-simplex. 
Other geometrical objects of interest here are the following: 

- two triangles tijh and tijk sharing the edge Cij identify an angle ^P[ij)hk which can be defined 
as TT minus the angle between the normals to the two triangles within the tetrahedron Tijhk ] 

- two tetrahedra Tijhk and Tijhi sharing the triangle tijh identify a dihedral angle 9(^ijh)ki which 
can be defined as tt minus the angle between the normals to the two tetrahedra within the 
4-simplex Y,ijhki- 
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The geometry of a 4-siinplex can be completely described giving the lengths Lij of its ten edges, 
provided that they satisfy a set of triangular inequalities and their higher dimensional analogues. 
Following [53], we introduce the {n + 2) x {n + 2) matrix 



^Ji.n+l ) — 
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The area of the triangle ^123 can be expressed in terms of the edge lengths L12, ii3, -^23 as 

1 



1123 



V- dot M2(Ll2, il3, L23) 



(118) 



an expression which corresponds to the well-known formula of Heron. Similarly, the volume of the 
tetrahedron T1234 can be written as 



V1234 — 



1 



3!2^ 



■\/ + det Mz{Li2, ii3, Li4, L23, L24, L34) 



(119) 



which corresponds to Tartaglia's formula. Finally, for the 4-volume of the 4-simplex we have 
1 



'12345 



4! 2^ 



-■\/— det M4(Li2, Li3, Li4, L15, L23, ^24, L25, ^34, ^35, -^45) 



(120) 



Triangular, tetrahedral and pentachoral inequalities on the edge lengths correspond to the require- 
ment that the areas of all the ten triangles, the volumes of the five tetrahedra and the 4-volume of 
the 4-simplex are given by positive real quantities when expressed in terms of the edge lengths using 
formulas (jllSp . (|119p and (|120p . Angles and dihedral angles can be expressed in terms of areas, 
volumes and the 4-volume. The angle (^(12)34 between the triangles ^123 and ii24 can be written in 
terms of the lengths of the six edges of the tetrahedron T1234 they belong to using the following 
formula 

3 -^12 

</5(12)34(il2, il3j -^14, -Z>23j -^24, -^^34) = Sin^ Vi234- (121) 



2 ^123^124 

Similarly, the dihedral angle ^(123)45 between the two tetrahedra ri234 and T1235 can be written in 
terms of the lengths of the ten edges of the 4-simplcx and is given by 



^(123)45 (-^u) — sin ^ - 



A 



123 



3 V1234V1235 



^1 



2345 



(122) 



As within a 4-simplex for each triangle tijh there is only one dihedral angle 0(^ijh)kii we use the 
simpler notation 0(ij7i) to identify it. A geometrical property the dihedral angles of a 4-simplex 
satisfy is the following 

Ahkl{Lmn) —^7 {Lrnn) — . (123) 



l<h<fe<;<5 



It goes under the name of Schlafli differential identity (see the appendix of ^ for a geometrical 
proof). 
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We call regular a 4-simplex which has all its 10 edge-lengths equal. Using Bang's theorem, Baez 
et al. [22] show that a non-degenerate 4-simplex which has all its 10 triangle-areas equal is neces- 
sarily regular. Now we consider fluctuations of edge-lengths SLmn around the regular configuration 
with lengths Lq and show that the nearly-regular 4-simplex they define can be described in terms 
of triangle-area fluctuations SAhki around the regular configuration with areas Aq — ^L^. In fact 
the Jacobian matrix U from length- fluctuations to area- fluctuations (see equation (|106p ) can be 
computed using formula ()118|) and turns out to be invertible. We call V its inverse. Choosing the 
following ordering for the fluctuation vectors 

{SL,nn} = {SLi2,SLi3, SL^, 5Li5,5L23, 5L24, 5L25, 5LzA, SLi^) , 

{SAijh} = ((5^345, M245, <^^235, (5^234, <5^145 , <^^135 , (5^134, SAi25,5Ai24, M123) , 

the matrices U and V turn out to be symmetric and have the following coefficient^^: 
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(124) 



(125) 



'^The matrices ||124| |. ||125| I and l|73|I . H108ll all have the same structure which corresponds to the following matrix 



/ a 



K{a,b, c) 



V 



b 
b 
c 
b 
b 
b 
b 

a J 



with specific coefFeicients a, b, c. Here we collect some useful formulae for the eigenvalues of a matrix of this form: 

Aq = a + db + 3c (non-deg.) 

Xl = a + b-2c (4-fold deg.) 

A2=a-26 + c (5-fold deg.) 
The eigenvalue Ao is associated to the eigenvector vq = (1, 1, 1, 1, 1, 1, 1, 1, 1, 1). 
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Using the matrix (|125p . the geometry of a nearly-regular 4-simplex can be completely described 
in terms of area- fluctuations: for instance the fluctuation of the volume of the tetrahedron T1234 is 
given by equation 

^^1234 = E ( E ^^(^0) ^i:mhki)) SAf,,i (126) 

l</i<fc</<5 l<i<j<i 

in terms of the ten area-fluctuations. 
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